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The fc-point correlation functions of the Gaussian Random Matrix Ensembles are certain deter- 
minants of functions which depend on only two arguments. They are referred to as kernels, since 
they are the building blocks of all correlations. We show that the kernels are obtained, for arbi- 
trary level number, directly from supermatrix models for one-point functions. More precisely, the 
generating functions of the one-point functions are equivalent to the kernels. This is surprising, 
because it implies that already the one-point generating function holds essential information about 
the fc-point correlations. This also establishes a link to the averaged ratios of spectral determinants, 
'. i.e. of characteristic polynomials. 

^ ' PACS numbers: 05.45.Mt, 03.65.Nk, 05.30.-d 
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■ I. INTRODUCTION 



Random Matrix Theory (RMT) allows one to model a rich variety of complex systemsiiSi^. The Gaussian Unitary 
Ensemble (GUE) of random matrices is used in the absence of time reversal invariance. The Gaussian Orthogonal 
Ensemble (GOE) and the Gaussian Symplectic Ensemble (GSE) apply if time reversal invariance holds and if the 
, levels are not or are Kramers degenerate, respectively. These three cases GOE, GUE and GSE are labeled by the 
OO Dyson index /3 = 1,2,4. Supersymmetry^ often yields a considerably clearer insight into the structure of random 
T-H , matrix models. This is so, because supersymmetry drastically reduces the number of degrees of freedom without 
' giving away any information contained in the model. Thus, supersymmetry removes a certain kind of redundancy. 
' Loosely speaking, the supersymmetric formulation plays for the random matrix model the role of an "irreducible 
. representation". 

^ In this contribution, we present an unexpected direct connection between the fc-point correlation functions of 
■~»>^ , the three Gaussian Ensembles and the generating functions of the one-point functions. The fc-point correlation 
functions are determinants (GUE) or quaternion determinants (GOE,GSE)ii. All entries of these determinants are 
II , fully specified by one function of two energy arguments. Because of their fundamental importance, these functions 
' are referred to as kernels. In the supersymmetric formulation, generating functions are used which, upon derivative 
I with respect to source variables, yield the fc-point correlation functions. Especially the generating function of the 
■ one-point correlation function, or rather one-point function, depends on an energy and on a somewhat unphysical 
' source variable. The source variable is needed to break a symmetry in the supersymmetric matrix model. Hence, the 
I> . number of dependend variables is the same for the generating function of the one-point function and for the kernels. 

We show in the sequel that, surprisingly, the generating function of the one-point function is fully equivalent to the 
^ , kernels. This is true for all three Gaussian Ensembles and for arbitrary level numbers. Thus, a fundamental link is 
' established between the one-point functions and the fc-point correlation functions. 

The article is organized as follows. For the convenience of the reader, we briefly compile the relevant formulae for 
the random matrix correlation functions and kernels in Section ^1 We follow closely Mehta's book"'^. In Section ITTTI 
we present our main results and discuss implications. The derivations are performed in Section llVI Summary and 
conclusion are given in Section El 



X 



II. RANDOM MATRIX CORRELATION FUNCTIONS AND KERNELS 

The fc-point correlation functions Rk{xi, . . . , Xk) are the probability densities to find fc energies at positions 
xi,. . . ,Xk, regardless of labehng. They can be written as averages over a probability density Pn{H) of a TV x TV 
Hamilton matrix i?, 

k 

Rk{xi,...,Xk)= PNiH)l[trSixp-H)d[H] , (2.1) 

p=i 

where d[H] is the volume element, that is, the product of the differentials of all independent variables in H. We 
mention in passing that this definition contains contributions proportional to S{xp — Xq). However, as this issue is not 
important here, we ignore it and refer the reader to the discussion of those details in Ref^. 
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For the Gaussian Ensembles, the correlation functions have a remarkable determinant structurei. All knowledge 

k 



needed to construct the full function R^^\xi, . . .,Xk) is contained in one single function, the kernel, which depends 



on two energy arguments. In the case of the GUE (/? = 2) for TV levels one has 

.,Xk)= det K^^\xp,Xq) 



where the kernel is given by 



N~l 



K)^'{Xp,Xq) = ^ tpniXp)tpniXq) . 



n=0 



Here, ipniz) denotes the oscillator wave function 



:exp ( -y ) Hn{z) 



(2.2) 



(2.3) 



(2.4) 



and Hn{z) is the Hermite polynomial of order . 

Due to the additional symmetries, the corresponding expressions for the GOE (/3 = 1) and the GSE {(3 — 4) are 
more involved. For the GOE, the kernel is given by 



(2.5) 



K^j^\xp,Xq) ^ K^jP{Xp,Xq) + ^ —IfiN-liXp) j e{Xq - z)ip N {z)dz + ON (Xp) , 



where K)^' {xp, Xq) is the GUE kernel and the function 



N odd 
N even 



(2.6) 



enters. We also use the notation 



e{z) = -sign (z) . 



In the case of the GSE, the kernel reads 



if (4). X _ —if 



2Xp , V2x 



2N+1 



ip2N{V2Xp) 



+ 00 



e{xq - z)ip2N+i{V2z)dz , 



(2.7) 



(2.8) 



the first two terms are the same as in the GOE kernel, but for 2N + 1 levels. However, the function aN{x) does not 
appear. It is convenient to scale the energy arguments with \/2- 

The ordinary determinant for the GUE correlation functions is replaced by quaternion determinants. One has for 
the GOE 



Ri'\xi, 



and, similarly, for the GSE 

Ri'\xu 



,Xk) 



qdet 



, Xfc) = qdet 



{Xp,Xq) DKj^ (Xp,Xq) 
JKj^ {Xp^Xq) -^J^ {Xq^Xp^ 



^ N (Xpj Xq) DKpj {Xpj Xq) 



I K^j^^ (^Xp, Xq) K'^^(^Xq,Xp) 



(4), 



(2.9) 



(2.10) 



Here, D,I and J are certain derivative and integral operators, respectively. We write Kj^ {xp,Xq) for the GUE 
kernel which Mehta denotes by K^ixp, Xq)^-. Mehta works with the kernel Sn{Xp, Xq) for the GOE and the GSE. This 



function S'Ar(xp, a;,) is our K^\xp, Xq) without the function ajqixp). We decided to introduce the kernels K)^' {xp, Xq) 

and K^^\xp, Xq), because the function aN^Xp) is only present in the GOE, but not in the GSE case. As we will show, 

the kernels K^\xp, Xq) are the ones that appear naturally in the supersymmetry context. More information on the 

relation between Mehta's kernels and the kernels K^\xp, Xq) and on how they enter the expressions (|2.9I) and (|2.10() 
for the correlation functions can be found in [XI 

In concluding this compilation, we underline once more that the knowledge of the three kernels suffices to build up 
all fc-point correlation functions for the three ensembles GOE, GUE and GSE. 



(1)/ 
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III. KERNELS, MATRIX INTEGRALS, GENERATING FUNCTIONS AND RANDOM MATRIX 

AVERAGES 



Surprisingly, one can obtain the kernels from the lowest dimensional one-point supermatrix models which reflect 
the appropriate symmetries. This and its implications state the main result of the present contribution. For the GUE, 
i.e. for /3 = 2, we have 



Kj^\xq, Xp) 



1 exp {xl/2 - xy2) 



where 



Im ( - / cxp (-trgcr^) detg~''(a - x-)d[a] - 1 ) , 



a X* 
A ib 



(3.1) 



(3.2) 



is a 2 X 2 Hermitean supermatrix. The entries a, b are real commuting and A is complex anticommuting. The 
energies are ordered in the diagonal matrix x = diag {xp, Xq). The variable Xp is supplemented with a small imaginary 
increment irj such that and X = diag {Xp , Xq). The corresponding result in the case of the GOE, i.e. for 

13=1, reads, for even and odd level number N , 



K%\:,^^^^) = if^^i^Zl-^ In, Q j cxp (-^trga^ ) detg-^/^(a - x-)d[a] - 1 
Finally, for the GSE, i.e. for /3 = 4, we have 



(Xq, Xp) — 



1 exp (, 



27r x„ 



i^lm j cxp (-trgCT^) dctg"^(CT - x-)d[a] - 1 



(3.3) 



(3.4) 



In the cases of the GOE and the GSE, cr is a 4 x 4 Hermitian supermatrix with an additional symmetry, often referred 
to as orthosymplectioii^. Explicitly, a reads 



y/ca y/cd X* 

\fcd y/cb /i* 

A /i \J~^w 

X* fi* 



-A 



cw 



(3.5) 



where c = 1 or c = —1 for GOE and GSE, respectively. Here the variables a,b,d and w are real commuting, while 
A, A* and /i, fi* are complex anticommuting. The diagonal matrix of the energy arguments now also has dimension 
4x4 and reads x = dia^g (xp, Xp, Xq, Xq). However, for brevity we always write x for the 2x2 and for the 4x4 
energy matrix. All formulae given here are exact for finite values of TV. Thus, as always, supersymmetry decouples 
the number of integrations to be done from the level number N. The number of integrations is fixed, while the level 
number N may take arbitrary values. The common form of the three results (|3.1|) to (|3.4|) is evident. The differences 
for the three Gaussian ensembles lie in the structure of the matrices a. Our definitions and notations are the standard 
onea^ and, in particular, they follow the definitions and notations of RefsiSifi. Thus, we may even formulate the 
results (|3.1|) to (|3.4|) in the compact form 



1 exp {-f {xl - xl) /2) 

77r Xp - Xq 



Im 



exp 



27 



detg-'^^/^l''l((7-x-)d[(7] - 1 



(3.6) 



where we introduced 7 = 1 for /3 = 1, 2 and 7 = — 2 for /3 = 4. 

In spite of its non-trivial character, the result (|3.1() is easily proven because it is an immediate consequence of 
an integral representation of the kernel K'"^\xp, Xq) which was found in Ref.^. We briefly sketch the derivation in 
Section Hvl In fact, expressions similar to Eq. (|3.1|l have already been used for a study involving the chiral GUEiS 
and for a certain generalization of the GUE^^. On the other hand, the proofs of the results (|3.3|l and (|3.4|l are more 
involved. They will also be given in Section llVI 

Formulae (|3.1|l to (|3.4|) are remarkable, because they establish a direct and previously unknown connection between 

the kernels and the generating functions z[^^ {x) of the one-point functions, 



R['\x^) 



d_ 

2|7|7r5Ji 



(3.7) 



Ji=0 
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We introduced the diagonal matrix x = diag {xi — Ji, + Ji) which contains the energy argument xi and the source 
variable Ji. To be consistent with the previous notation, we use Xi to denote the argument of the one-point functions. 
The one-point function is written as r[^\xi) = r[^\xi) + iR[^\xi) such that the level density is the imaginary 

part, Imi?^'^''(a;i) = Rf \xi). As is well known, the matrix integrals in the expressions H3.1|l to H3.4|l are precisely the 
generating functions, 



Zn^) = ^ / exp ( -T^trga^ ) detg-^^' ^^^-^^a - x-)d[a] (3.8) 



8^4 J --^ y 2|7 
with the appropriate supermatrices a. Thus, we arrive at 

1 ..^.2/0 .2/o^T^^^(^)-^^(o) 



K^j^\xq,Xp) = -exp(xp/2-x^/2)lm- 

(1), . 1 / 9 , 9 , N Z^^^\x) - zf^(O) 
K)^\x„xp) = -exp(a;2/2-x2/2)lm- ^ ^ ^ ^ ^ ^ 



K^j^\xg,Xp) = —exp{xl-xl)lm ' ^ ^ ^ ■ (3.9) 

For the GSE kernel K^\xq, Xp) the arguments of the exponential are interchanged with respect to the GUE kernel 
K^^\xq, Xp) and the GOE kernel K^^\xq,Xp). We notice that z'{^\x) depends on the two energies Xp and Xq. There 
is no source variable here. Moreover, we have Z^^^\q) — 1 due to the definition of the generating function. Again, we 
can write 



^ »p (2 (4 - .J)) i„, (3^0) 

which combines the three results (|3.9|l in a compact form. 

Formulae 1)3. 9|l state a close connection between the kernels and the generating functions. The kernels can be viewed 
as difference quotients of the generating functions at the two points x and 0. The crucial quantity is the difference 
Xp — Xq. By construction, the generating functions are unity whenever the two arguments degenerate. Thus, Z^^\x) 
moves away from unity as function of Xp — Xq. If one takes the limit Xq — > Xp, the difference quotient becomes the 
differential quotient (|3.7|l . This yields the kernels at Xp = Xq, that is, the level densities as function of the single 
remaining variable. In ^ we discuss extensions of the previous results if real parts contribute to the correlation 
functions. 

To further clarify the meaning of these findings, we rewrite the generating functions as averages over the original 
random matrices in ordinary space, 

with normalization constants C^v/j. Here, the matrices H parameterize the GOE, GUE and GSE for /3 = 1,2,4, 
respectively, as defined in Mehta's book^. Combining Eqs. (|3.9|) and H3.11|l . we see that the kernels themselves are, 
apart from factors, averages over the Gaussian ensembles. This is a surprising insight. According to the definition (|2.1ll . 
the A; -point correlation function of a Gaussian ensemble is one single matrix integral for a fixed value of k. The results 
presented here imply that this single average breaks up into products of averages. This is intimately related to the 
determinant structure, but it is a stronger statement because it identifies the determinant structure as stemming from 
the break up of the random matrix average. 

Furthermore, it follows from Eqs. 1)3. 9|l and 1)3.11(1 that the random matrix kernels are essentially an average over 
a ratio of spectral determinants, taken at the two different energies. This relates our findings to the presently much 
discussed issue of characteristic polynomials, spectral determinants and their moments, see Refs^Sii^ and references 
therein. For matrix dimension N — 2, the connection between averages over the ratio of two characteristic polynomials 
and the kernel was recently observed in Refii^ in the GOE case (/3 = 1). 



IV. KERNELS EXPRESSED AS EIGENVALUE INTEGRALS IN SUPERSPACES 



We prove the results in the previous Section by explicit calculation. Alternatively, one could try to employ Dyson's 
Brownian motion^ and its supersymmetric extensioiJi for the stationary case, i.e. for the pure ensembles. However, 
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this would still leave one with the problem of fixing the boundary conditions in an unambiguous way. Another strategy 
could consist in showing that the supermatrix models satisfy the same equations that the kernels obey, such as the 
convolution condition. Once more, one is confronted with some ambiguity. Thus, we believe that the most direct 
proof is probably an explicit caculation, but we certainly do not exclude that other direct proofs also exist. 

There are two possibilities to proceed with an explicit calculation. First, due to the small dimensions of the 
supermatrices in Eq. H4.1(l . one can expand the superdeterminants in the supermatrix models and integrate out the 
Grassmann variables by "brute force" . The resulting expressions are rather complicated and the calculations to 
follow are quite cumbersome. Second, one can introduce eigenvalue-angle coordinates and integrate in a first step 
over the supergroups and in a second one over the eigenvalues. We present this approach in the sequel because the 
eigenvalue integrals to be solved here are of a general type which will always appear in exact calculations involving 
supersymmetry. In particular, they will show up in generalizations of the present supermatrix models. Thus, we want 
to develop techniques for how to handle them. 

We denote the right hand sides of the formulae l|3.1() to H3.4|l by 

where, for notational convenience, the exponential functions are not included. For the same reason, we also split off 
a factor The strategy for the ensuing calculation is now the same in all three cases (3 = 1,2, 4. The supermatrix 
integral yields two contributions, an Efetov-Wegner termi& and an eigenvalue integral. The Efetov-Wegner term 
gives a constant, canceling the unity to be subtracted in the bracket. To obtain the eigenvalue integral, we shift 
the integration matrix a by the diagonal matrix x, introduce eigenvalue-angle coordinates a = u~^su for the shifted 
integration matrix and perform the group integral over the diagonalizing matrix u. For (3 = 2, the group integral is 
the well known supersymmetric Itzykson-Zuber integral^. It is trivial here, because u is only a 2 x 2 supermatrix. For 
/3 = 1,4 the group integral is over 4x4 supermatrices u and thus non-trivial. In both cases /3 = 1,4, it is the same 
group integral which can be viewed as a supersymmetric extension of Gelfand's spherical functions. This integral was 
first calculated in Refsn^. (Responding to some confusion in the literature, we mention in passing that the integrals 
involving 8x8 supermatrices were also first solved in these references. However, in the present case, the 4x4 case 
suffices.) After the group integrations, we are left with integrals over the eigenvalues s = diag(si,?S2) for (3 = 2, 
s = diag (sii, S21, JS2, is2) foi (3 = 1 and s = diag (si, si, isi2, JS22) for /3 = 4. These eigenvalue integrals are solved 
in the sequel. We notice that the eigenvalues in the fermion-fermion block of the supermatrix a carry an imaginary 
unit. It is due to a Wick rotation which is necessary to ensure convergence of the supermatrix integrals. Thus, this 
imaginary unit has to be ignored when taking the imaginary part on the right hand side of Eq. (|4.1fl . 



A. Gaussian Unitary Ensemble 



The energy difference Xp — Xq drops out and the matrix integral reduces to the double integral 

L^^\xp,Xq) = -^ r r '^''^'^ exp (-(si + xpf + {iS2 + xqf) {iS2flxn . (4.2) 

This coincides with the double integral found in Ref:-. In this reference, it was denoted by CN{xp,Xq). However, 
it is important to notice that the double integral resulted in Rcf.^ from calculating the /c-point correlation function 
for arbitrary k, that is, from a group integral over a 2k x 2k unitary supermatrix. In Ref.^, the double integral was 
already evaluated and it was shown that 

K^N^ (xp, Xq) = exp (4/2 - xl/2) L^^ (xp, Xq) . (4.3) 

This proves formula 1)3. l|l . 



B. Gaussian Orthogonal Ensemble 

The orthogonal case has a much more complicated structure. From Refsi2i^, a triple integral results, 

^(1), ^ _ f^"^ [^°° [^°^ \sil - S2l\dsilds2ldS2 

\-^p' ■^q) 



S'T^^ J-oc. J-oo J-oo {Sll- iS2Y{s21- iS2Y 
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exp + XpY - ^(*21 + Xpf + {iS2 + Xqf 

{2{Xp - Xq){sii - iS2){s21 - 182) + (sn + S21 - 2iS2)) 

Again, the energy difference the denominator has been canceled. The technical difficulty to overcome is 

twofold. First, the integration variables are coupled in an inconvenient way, even involving an absolute value. This 
happens in the term in front of the Gaussians which stems from the Jacobian or Berezinian of the eigenvalue-angle 
coordinates. Second, the last term contains an imaginary part of the product of two singularities. It cannot be 
interpreted as the product of imaginary parts which would simply yield derivatives of 5 functions. Both difficulties 
can be circumvented by observing that the GOE kernel contains, according to Eq. (|2.5|) . the GUE kernel. Thus, we 
split off the latter. To this end, we reformulate the GUE kernel H4.2|l as a triple integral. We write si = sn and 
introduce a dummy integration over the variable S21 by multiplying the integrand with the function i5(sii — S21). 

A2), ^ J_ /■+°° Sjsn - S2i)dsnds2ids2 

\Xpi Xq) 



2tT^ J-oo J-oo J-oo (Sll - «S2)(S21 - iS2) 

exp ^-^(■511 + ^p)"^ - ^(^21 + Xp)'^ + {iS2 + Xqf 

(-" + --^'^--)(--)"^'" (.rj^42ir/- (4.5) 



In order to subtract this formula from Eq. (|4.4|l we have to do some integrations by parts in both expressions. More 
precisely, we use 

5(.n-.S2i) = i(^-^)e(.n-.2i) (4.6) 

in Eq. 14.5|l . This procedure casts the integrand in Eq. H4.5|) into the adequate form to be subtracted from the left 
hand side of Eq. (|4.4|) . It is also convenient to do an integration by parts in Eq. 14. 4|) using 



n / \ 9 d d ^ 

2 (Xq - Xp) - — h 2lS2 - sii - S21 

diS2 osii ds2i 



exp ( -^(sii + Xpf - i(s2i + Xpf + [is2 + Xqf ) = . (4.7) 



With these adjustments we can subtract l''^ (xp,Xq) from L^^\xp,Xq) and obtain 



iXp,Xq) — Lj^ (Xp,Xq) L (Xp^Xq) 



— / / - S2l\dsiidS21 

^ J — 00 J —00 



N 
8^ 



exp ( --^(sii + Xpf - i(s2i + Xpf ) Im ■ 



2' 2' J (sn)^/^+Hs2i)^/^+^ 

+00 

ds2exp{is2 + Xq) {is2) ~ . (4.8) 

-00 

This expression decouples the three-dimensional integral into a product of a two-dimensional integral and a one- 
dimensional integral. Furthermore, the left hand side of Eq. (|4.8|) factorizes into a product of functions each depending 
only on one energy argument. Therefore, we can write 

M^^\xp,Xq) = ^J^\xp)^^^\xq) (4.9) 

The function i/;^'' is simply an integral representation for the Hermite polynomial. The integration of requires 
more effort. The relations 

^n\^p) - (y + 1) '^^Ul^^p) = 4^ [~N + iy. exp(-4) (4.10) 
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are used. The second formula above was derived by using 

^exp(.^)lnry_^^ 



HM-p) = '-^^ exp (4) Im / ^^^V_W+r' ' '^^ (4-11) 



and by the introduction of a dummy variable similar to Eq. I|4.5|) . By combining these relations, we obtain 
^^^(^^p) = -exp(-xg/2) ( ^^^"^^ y e(xp-t)i/^(t)exp(^-^^di + c(^^^ 

= ^^^^^i^HA.-iK). (4.12) 
The form of the above equations are the expected ones, and the remaining problem is to calculate the integration 



constant c^-*. This tedious calculation is performed in [21 We find 



(1) _ J if iV even 



"^JV - ^ _47r2^/ViV!! if iV odd " ^^'^^^ 



This non-vanishing constant gives rise to a contribution which is identified with the function cjat defined in Eq. 

Now, Eq. H4.12|l is rewritten in terms of the oscillator wave functions ipn defined in H2.4|) and compared to H2.3|l and 
()2.5f) . We then have 

M'^\xp,Xg)^exp{xy2-xl/2) (^K'}^\x,,Xp) - K^j^\x„Xp)) (4.14) 
for all values of N. This proves Eq. H3.3|) . 



C. Gaussian Symplectic Ensemble 

For the GSE the structure of the supermatrix a is almost the same as for the GOE. The group integral found in 
Refs.^'^ can be applied again. However, boson-boson block and fermion-fermion block are interchanged with respect 
to the GOE. As a consequence the imaginary unit now comes in front of the integration variables sii,S2i and the 
contribution of the superdeterminant in Eq. (|4.4|) is inverted. With an additional rescaling a — > cr/V2 we obtain 

+ 00 /' + 00 /' + 00 



(4) 1 r-l-oo /• 



\sil - S2l\dsiids2ldS2 
2 



(iSll - S2Y{iS2l ~ S2Y 

exp f - {isii + \/2xj^ + - {is2i + V^Xj^ - + V^x^ j 

^2 (^\f2xq ~ \/2xj^ (isii - S2)(is2i - S2) + (isii + is2i - 2s2)) («Siiis2i)^ Im ^^^^^n ' (4-15) 
Now we can apply the same method as in the case of the GOE. We arrive at the decomposition 

L''^\xp,xq) = 4w(V2x„V2xp) +a4^^(xp,x,) (4.16) 
M^^\xp,x,) = '^ujf{xp),l,f{x,). (4.17) 

The functions oj^'' and V'^'' are now given by 



ds2 exp ^- \ s2 + \/2a;g j j Im 



The function il}'^ (xq) is easily evaluated by making use of the identity Eq. H4.11|l for Hermite polynomials 



y dsiids2i|sii - S21I exp f - ^isii + y2xpj + - (^«S2i + V2xp j j («siiis2i)^"^ (4.18) 



V'^^(^,)= (^exp(-2a;2)i/2jv(%/2x,) . (4.19) 
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The integration of Lu^j^\xp) is a little more tricky. The procedure follows ideas analogous to those used for the 

calculation of ^^•'(xp) for the GOE in Sec. (|IV B|l andlCl An integration constant occurs in this case as well. It can 
be fixed in the same manner as for the GOE. Here, however, it vanishes for all values of N. One finds 

^n\^p) = -^^^^p{4) J e(\/2xp-t)exp(-iV2)i?2W-i(i)dt . (4.20) 
Inserting these results into Eq. H4.16|l and expressing everything in terms of oscillator wave functions one arrives at 

C2W-1 
Tl = 

\j—<P2N{^Xq) j e[y2xp~t)<p2N-i{t)dt\ . (4.21) 
This is almost the final result. We use the integration formula 



2N r°° , . _ , _ /2iV + 1 



+ 00 



e{x -t)ip2N-\it)(ti. = •f2NWl + \\ / £{x-t)ip2N+\[t)dt (4.22) 



to obtain 



L^^N^^v^^q) = exp(a;p-a;^) ^i^^^^^i (^\/2xp,y2x,^ 

+ Y ^ ^2N{V2xq) J e[V2xp-tjip2N+iit)dt] . (4.23) 

This is exactly our assertion (|3.4|) . 



V. SUMMARY AND CONCLUSION 



We showed that the generating functions for the one-point functions yield directly the kernels of the correlation 
functions in RMT. This is tantamount to saying that the kernels are given by the lowest dimensional supermatrix 
models. We proved this by explicit calculations for the Gaussian Ensembles GOE, GUE and GSE. Recent results for 
supergroup integrals enter our derivation. We develop new techniques for integrals over eigenvalues of supermatrices. 
This was another reason for us to prove our results by explicit calculation. 

The equivalence between kernels and the generating functions for the one-point functions is an unexpected, surpris- 
ing insight. The generating functions contain a symmetry breaking, the source term. Our results demonstrate that 
this symmetry breaking is intimately related to the correlations themselves. Hence, the generating functions comprise 
much more information than just the one-point functions. Technically, this becomes apparent in the fact that the 
source variable adds, with different signs, to the energy variable. There are effectively two energy arguments which 
are then identified with those of the kernels. The kernels are obtained as the difference quotient of the generating 
functions. In the limit of the differential quotient, one finds the well known relation between one-point functions and 
their generating functions. 

All our results hold for arbitrary level number. Among other things, this opens yet another possibility to calculate 
the kernels in the limit of large level number on the local scale. Here, one can do that by a saddlepoint approximation 
of the one-point supermatrix models. No Goldstone modes are present and one finds the kernels for all correlations 
from the saddlepoints. As this is a straightforward exercise, we have not presented it in this contribution. 

Our findings are likely to have further extensions. From ReffJA, one easily concludes that the structural relation 
we observed carries over, for /3 = 2, to models in which a fixed matrix is added to the random matrices. Further 
investigations are in progress for the cases (3=1 and /3 = 4. Our results could have relevance for field theory as well. 
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APPENDIX A: MEHTA'S KERNELS AND THE KERNELS K'-^\xp,Xq) 

The fundamental piece in Mehta's notation is the function 



Therefore in his expressions for the GOE correlation functions a n (x) appears as an independent quantity 



+00 



R'j^\xi, . .. ,Xk) ^ qdet 



SN{Xp,Xq) + a{Xp) DSN{Xp,Xq) 

JSN{Xp,Xq) SNiXq,Xp) + a{Xq) 



>,q—l,...,k 



In the GSE result, the function a{z) does not appear, 



rI^\xi, ...,Xk) = qdet 



1 
71 



S2N+liV^Xp,V2Xq) DS2N+liV2Xp,V^Xq) 
IS2N+l{V^Xp,\/2Xq) S2N+l{V2Xq,\/2Xp) 



(Al) 



(A2) 



(A3) 



p,q—l,...^k 



The operators D,I and J are defined as acting on the function S{xp,Xq) only, 

DS]\[{Xp,Xq) = —-J^Sn{Xp, Xq) 



ISNiXp,Xq) = j dte{Xp - t)SN{t,Xq) 

JSN{xp,Xq) — ISN{xp,Xq)+ / a{t)dt — j a{t)dt + e{xp — Xq) 



(A4) 

Jo Jo 

In our approach, the kernels k''^\ i.e. the complete upper left entries of the 2x2 matrices in Eqs. (|A2I) and HA3|I . are 
the fundamental quantities, rather than Mehta's kernels. Therefore the operators in the off diagonal elements should 
also be defined as acting on K^'^^K This is accomplished by the following definitions 



1 

2 V dx,, 
1 
2 



-KN{Xq,Xp) - {Xp ^ Xq) 



dte{Xp -t)K^j^\t,Xq) - (Xp ^ Xq) 



(A5) 



JK}^'{Xp,Xq) = IKNiXp,Xq) + e{Xp - Xq) . 

With these definitions our expressions for the correlation functions, Eqs. 12.9|l and H2.10|l . are identical with Eqs. HA2|I 
and ljA3|l . This is easily verified by using that DS^ixp, Xq) and ISN{xp,Xq) are antisymmetric in their arguments^. 

We just remark that the simplicity of the definitions (|A5|I is another strong hint that the functions K^j^\xp, Xq) rather 
than S{xp,Xq) are the fundamental quantities. 



APPENDIX B: REAL PART CONTRIBUTIONS TO THE CORRELATION FUNCTIONS 

The correlation functions Rk{xi, . . . ,Xk) in classical RMT are, according to Eq. (|2.1|l . averages involving only the 
imaginary parts of the Green functions. Including the real parts, one has the more general correlation functions 

i?fe(xi,...,Xfc) = 4 / PNiH)l[tT—^d[H] . (Bl) 

We use the notation of Ref."'^^, cf. Eq. (|3.7|) . As for the definition (|2.1|) . we ignore contributions proportional to 
S{xp — Xq). In the case of the GUE, it has been shown in Ref.^ that the functions (jBl|) also have a determinant 
structure. We conjecture that the quaternion determinant structure carries over to the GOE and the GSE cases, too. 
The corresponding kernels K^\xq, Xp) are generalizations of the kernels K^^\xq, Xp). We expect that they are given 
by 

1 e^v{l{xl-xl)/2) ( (3^ /3 ,,„2\,,„-W2|7l, 



K^^>{xq^xp) = - ' (^^yexp(^-^trga^)detg-^"/-'^'(a-x-Ha]-l) , (B2) 
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such that Eq. (|3.6|l results when taking the imaginary part. In the GUE case /3 = 2, formulae ljB2p is a immediate 
consequence of Ref.^. For the GOE and GSE cases (3 = 1 and /3 = 4, formulae HB2|I states a conjecture. The kernel 

^ (2) (2) ^ 

{xq, Xp) follows from Kj^ {xq, Xp) by simply replacing one of the oscillator wave functions ipn{z) with (^„(z). The 
latter function combines the two independent solutions of the oscillator wave equation, i.e. the function ipn{z) and its 
Cauchy or Stiltjes transform. To the best of our knowledge, the relevance of those second solutions in an RMT context 
was first observed in Refiii. Again, we conjecture that these features also carry over to the kernels K^^\xq,Xp) and 
Kl^\xq, Xp) for /3 = 1 and /3 = 4. 



APPENDIX C: CALCULATION OF SOME INTEGRATION CONSTANTS 



Considering Eq. H4.10|l at Xp = 0, we obtain the recursion formula 
where 

bN= e{t)exp{^t^/2)HN{t)dt . (C2) 



The right hand side of Eq. (|(]l(l turns out to be zero for all N. This is easily seen for even N. For odd N, one has to 
employ Eq. H4.6|l and to integrate by parts. The recursion formula obtained in this way is equivalent to a result given 
by Mehta^ The remaining task is to find Cq^^ and c^^^ as starting values for an induction. 

We employ the implicit definition of uj[^\xp) according to Eqs. H4.8|l and H4.9f) . The difficulty is due to the 
singularities. For iV = 0, it suffices to use 



1 1 



(C3) 



■5ll^21 ■^11 '"'21 \^21 ^1 

A straightforward calculation and comparison with Eq. H4.12|l gives 

= 0. (C4) 

For = 1, the singular terms involve fractional exponents and the steps needed are more complicated. One can 
employ an integral representation of the F function, valid for arbitrary k. It yields 



1 



poo 

/ dtt'"'^ exp(-iis~^) 
Jo 



(C5) 



which moves the singularities into the exponent and decouples them from the power k. Moreover, we introduce new 
integration variables 

= and = ^— ^ with r = Tcosi9. (C6) 



All this leads to 



I Si 1 So 



' poo poo 

Jo Jo 



' ' ^°°drT^'=-iexp(-f (sri + s,-,)) 



sin^'^-i ^ exp '''^ cos ^ 



2(2fe-i)/2r(fc) 



^ drT2fc-iexp(^-y (s-^+s^-) 

'^(2fc-l)/2(T(sil - S2l)/2) 

(T(sn-S2i)/2)(2fe-i)/2 ■ 



(C7) 
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Here, we did the angular integral using the representation 

£ exp i^z cos z9) sin^~2 = 2(''-2)/2^p '^^^Hy^^ , (C8) 

for the Bessel function in d dimensions. We now insert Eq. (|C7|I into the implicit definition of uj[^\xp) according to 
Eqs. (|4.8() and (|4.9() . We also rotate the eigenvalues 

sii + S21 Sii - S21 

u = and V — (C9) 

and find, at Xp = 0, 

u;«(0) = Im I^Jdsn d.^i^^-^ exp (^-1 + 4)) 

= Im— / rfu / / dTTe{v)exp{-u'^ - v'^) cxp{-iTu-)Ji{Tv) 

r(3/2) Jq 

= -87r + 4\/27r (CIO) 
with the step function e{v) defined in Eq. p. 7(1 . From Eq. H4.12|l . we also have 

wW(0) = -87r^ exp^-it^j^di-cj^^ =-87r-4^^ (CU) 

Hence, combining the last two formulae, we obtain 

c[^'^ = -AV2tt (C12) 

for the case = 1. 

Thus, we can now use the recursion (|C1|I and finally arrive at the result (|4.13|) . 



^ M.L. Mehta, Random Matrices, 2nd edition, (Academic Press, San Diego 1990). 
^ F. Haake, Quantum Signatures of Chaos, 2nd edition, (Springer Verlag, Berlin 2001). 
3 T. Guhr, A. Miiller-Groeling and H.A. Weidenmiiller, Phys. Rep. 299, 189 (1998). 
* K. Efetov, Adv. Phys. 32, 53 (1983). 

^ J.J.M. Verbaarschot, H.A. Weidenmiiller and M.R. Zirnbauer, Phys. Rep. 129, 367 (1985). 

® M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, (Dover Publications New York 1970). 

T. Guhr and H. Kohler, math-ph/0012047 
8 T. Guhr and H. Kohler, J. Math. Phys. 43, 2741 (2002) 
® T. Guhr, J. Math. Phys. 32, 336 (1991). 
^° T. Guhr and T. Wettig, Nucl. Phys. B506, 589 (1997). 

T. Guhr and H.J. Stockmann, Journal of Physics A, in press 
E. Strahov and Y. Fyodorov, math-ph/0210010 
J. Baik, P. Deift and E. Strahov, math-ph/0304016 
^'^ E. Brezin and S. Hikami, J. Phys. A36, 711 (2003). 

T. Guhr, Ann. Phys. (NY) 250, 145 (1996). 
^® M.J. Rothstein, Trans. Am. Math. Soc. 299, 387 (1987). 
T. Guhr, Nucl. Phys. A560, 223 (1993). 



